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ABSTRACT 

CCOs are X— ray sources lying close the center of supernova remnants, with inferred 
values of the surface magnetic fields significantly lower (< 10 11 G) than those of 
standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting 
the first 2D simulations of the submergence and reemergence of the magnetic field in 
the crust of a neutron star. A post-supernova accretion stage of about 1O _4 -1O _3 M 
over a vast region of the surface is required to bury the magnetic field into the inner 
crust. When accretion stops, the field reemerges on a typical timescale of 1-100 kyr, 
depending on the submergence conditions. After this stage, the surface magnetic field 
is restored close to its birth values. A possible observable consequence of the hidden 
magnetic field is the anisotropy of the surface temperature distribution, in agreement 
with observations of several of these sources. We conclude that the hidden magnetic 
field model is viable as alternative to the anti-magnetar scenario, and it could provide 
the missing link between CCOs and the other classes of isolated neutron stars. 
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1 INTRODUCTION 
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The handful of reported Central Compact Objects ( CCOs ) 
forms a class of X-ray sources (see Ide Lucal ([2008); 
iHalpern fe Gotthelj (2010) for recent reviews) , located close 
to the center of ^kyr old supernova remnants. CCOs are 
supposed to be young, isolated, radio-quiet neutron stars 
(NSs). They show very stable, thermal- like spectra, with 
hints of temperature anisotropics, in terms of large pulsed 
fraction, or small emitting regions for hot (0.2-0.4 keV) 
blackbody components. 

The period is known for only three cases: P — 
424 ms for IE 1207.4-5209 in G296. 5+10.0 (|Zavlin et al.1 
l200Ch. P 112 ms for RX J0822.0-4300 in Puppis A 
(|Gotthelf fe Halpernl 120091 ). a nd P 105 ms for CXOU 
J185238.6+004020 in Kes 79 (jGotthelf et al.1 l2005h (here- 
aft er IE 1207, Puppis A and K es 79, respectively). For Kes 
79, IHalpern fe Gotthelj (|2010h reported a per iod derivative 
of P ~ 8.7 x 10" 18 ss _1 ; for IE 1207, IHalpern fe Gotthelj 
(2011) give two equally good timing solutions, with P — 
2.23 x 10" 17 ss" 1 and P = 1.27 x 10" 16 ss" 1 . Only an up- 
per limit P < 3.5 x IP" 1 * 3 ss" 1 is available for the CCO in 
Puppis A jGotthelf et alJl20ld : lde Luca et alJl2012h . Apply- 
ing the classical dipole-braking formula gives an estimate for 
the dipolar component of the external magnetic field (MF) 
of B p ~ 10 10 - 10 11 G (at the pole). This low value has 
led to the interpretation of CCOs as "ant i- magnet ars" , i.e., 
NSs born with very low MFs, which have not been ampli- 
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fied by dynamo actio n due to their slow rotation at birth 
(|Bonanno et alJl2005h . 

An alternative scenario to the "anti-magnetar" model 
that explains the low values of B p is to consider the fallback 
of the debris of the supernova explosion onto the newborn 
NS. During a time interval of few months after the explosion, 
it accretes material from the reversed shock, at a rate far 



superior than Eddington limit (IColgatel 1971 1 Blondinl 1986; 
IChevalieil[l989l : lHouck fe Chevaliertil99ll : ICofpi et al.lll996l ). 
This episode of hypercritical accretion could bury the MF 
into the NS crust, resulting in an external MF (responsible 
for the spindown of the star) much lower than the internal 
"hidden" MF. 

When accretion stops, the screening currents localized 
in the outer crust are dissipated on Ohmic timescales and the 
MF eventually reemerges. The proc ess of reemergence has 
been e xplored in past pioneer works ([Young &; Chanmugaml 
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simplified ID models and always restricted to dipolar fields. 
It was found that, depending on the depth at which the MF 
is submerged (the submergence depth), it diffuses back to 
the surface on radically different timescales 10 3 — 10 10 yr. 
Thus, the key issue is to understand the submergence pro- 
cess and how deep can one expect to push the field during 
the accretion stage. This latter issue has only been studie d 
(also in ID and for dipolar fields) bv lGeppert et all ([l999), 
in the context of SN 1987A. They conclude that the sub- 
mergence depth depends essen tially on the total accreted 
mass. More recently, IHoI (2011) has revisited the same sce- 
nario in the context of CCOs, using a ID cooling code and 
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Figure 1. Initial configuration of model A (left) and model B 
(right). Poloidal field lines (solid) smoothly match a vacuum so- 
lution outside the star. Upper panels: global configuration, the 
dashed line indicates the core-crust interface (denoted by R c )\ 
the color scale represents the toroidal MF intensity. Lower pan- 
els: planar representation of the MF in the crust (r E [R C ,R]) : 
as a function of the polar angle 0; the grey scale is normalized to 
the maximum value of Bt, and the white/black regions indicate 
positive/negative values. 

studying the reemergence phase of a buried, purely dipo- 
lar field, with similar conclusions to previous works. The 
hidden magnetic field scenario has also been proposed by 
IShabaltas fcTal (|2012h for the CCO in Kes 79. They can 
explain the observed high pulsed fraction (f p « 60%) of the 
X— ray emission with a sub-surface MF of ~ 10 14 G, that 
causes the required temperature anisotropy. 

In this paper, we further explore the viability of this 
scenario with improved calculations that can account for 
temperature anisotropics. We present the first results from 
2D simulations of the submerge nce and rediffusion of the 
MF, using the recent extension (IVigano et al.l|2012r) of the 
magneto-thermal evolution code of Pons et al. (2009b- The 
new code is able to follow the coupled long-term evolution of 
MF and temperature including the non-linear Hall term in 
the induction equation. This allows to follow the evolution 
for any MF geometry (not restricted to dipoles) and includes 
state-of-the-art microphysical inputs. 



2 BASIC EQUATIONS AND INITIAL MODEL 

We assume spherical symmetry for our background star 
(very small deformations induced by the MF or rotation 
are neglected) with redshift and other metric factors only 
depending on the radial coordinate, and we adopt the stan- 
dard static metric 

ds 2 = _ c 2 e ^(r) dt 2 + e 2X(r) dr 2 + ^2 > (1) 

where e~ x ^ = (1 — 2Gra(r)/c 2 r) 1//2 . The density, pressure, 
enclosed mass m(r), and v(r) profiles are obtained by solv- 



ing the hydrostatic equilibri um equations. We consider a 
1. 4 My NS (see S ection 2 in IPons et ail (|2009h . Section 4 
of Aguilera et al. (2008) and references within for all de- 
tails about the equation of state and other microphysical in- 
put), with two different initial configuration, shown in Fig.[T] 
Model A has an initial dipolar component of the MF ex- 
tended to the core, with a toroidal co mponent to mimic 
a typical MHD equili brium configuration (Ciolf i et al.l 2009; 
lLander fe Jon es 2009), in which the toroidal field is confined 
in the internal region denned by the last closed poloidal field 
line. Model B is a crustal confined magnetic field, with an 
extended (in the angular direction), quadrupolar toroidal 
field. The upper panels of Fig. [T]show the global initial MF 
geometry, which in both cases matches with an external vac- 
uum dipolar solution, with B® — 10 14 G. The toroidal field 
is Bt = 10 14 G at its maximum. As we are interested in the 
evolution in the crustal region, for visualisation purposes we 
show in the lower panels (and hereafter) a planar represen- 
tation of the MF in the crust. 

When the te mperature of the NS drop s below the melt- 
ing temperature (|Baiko &; Yakovlevl Il996h a solid crust is 
formed. Within few days, the inner crust up to the neutron 
drip density is crystallized, while the outer crust is formed 
within weeks or mo nths (|van Riper! Il99ll ; lAguilera et al.l 
120081 : iHo et al J 12013 ). The evolution of the MF in the solid 
crust is described by the Hall induction equation that, in- 
cluding relativist ic corrections, has the following form 

^ = -Vx[r|Vx (e"B) - v e x (e u §)] , (2) 

where 77 is the magnetic diffusivity, and v e is the "electron 
fluid" velocity given by 

i-f^Vx (e"B) (3) 
47ren e 

with n e being the electron number density. The differen- 
tial operator V associated to the spatial metric includes the 
corresponding metric factors, i.e. , e~ x ^ In the core, the 
field is practically frozen because of the high conductivity of 
nuclear matter. In the solid crust, in the absence of accre- 
tion, the evolution of magnetic field is governed by Ohmic 
dissipation (the first term on the right hand side of Eq. ([2])) 
and the Hall term (term proportional to v e ). 

To compare the relative importance of the two terms, 
we denote by ujBT e = cB / Anen e rj the so-called "magnetiza- 
tion parameter", usually found in the literature of plasma 
physics, where ujb is the gyro- frequency and r e the elec- 
tron relaxation time. Note that this quantity gives an es- 
timate of the ratio between the characteristic Ohmic dis- 
sipation timescale Td = L 2 /r] and Hall timescale th = 
A.7ien e L 2 / cB — Td / (w BT e ) • This ratio strongly depends on 
temperature. Typical values in neutron stars vary in a 
large range T d ~ 10 4 - 10 9 yr and T h ~ (10 4 - 10 6 )^^ 
yr. In previous works ([Pons &; Geppertl 120071 ; I Pons et all 
2009; Vigano et al. 20JJ)7 we have discussed in detail the 
timescales, the mathematical character Eq. ([2]) and the nu- 
merical metho d to solve it. 

Following (jGeppert et al.lll999h . we introduce a new ve- 
locity term (v e — >> v e +v a ) in the induction equation to model 
the effect of accretion. This new term represents the sinking 
velocity of the accreted matter as it piles up on top of the 
crust, with 
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Figure 2. Radial profile at the equator of the Ohmic timescale 
(solid) and Hall timescale (dashed) for initial model A, with B® = 
10 14 G and L = 1 km. 



where rh is the mass accretion rate per unit normal area, 
which in general depends on the latitude (0), and r is the 
unit vector in the radial direction. This expression is a result 
of the continuity equation (conservation of mass) together 
with the assumption that the accreted mass steadily piles 
up pushing the crust toward the interior without spreading 
(pure radial displacements, thus with rh constant at any 
depth within a vertical column). 

3 MAGNETO-ROTATIONAL EVOLUTION 

Our code follows the evolution of the MF inside the star. We 
impose as a boundary condition at the surfac e the continu- 
ous m atching with general vacuum solutions ([Vigano et al.l 
l2012h . From the values of the MF at the surface, we can 
reconstruct the external configuration. In particular, we can 
follow the evolution of the external dipolar component, B p , 
responsible for the spin-down of the star. 

3.1 Submergence of the magnetic field. 

We study the accretion phase by exploring the sensitivity of 
the results to two key parameters: the total accreted mass, 
M a , and the angular dependence of rh(0) in Eq. J4}. Neither 
the final geometry nor the submergence depth at the end of 
the accretion stage depend on the duration and time depen- 
dence of the fallback accretion period that follows a super- 
nova explosion, because of its short duration (~ months), 
much shorter than the other relevant timescales shown in 
Fig. El 

It should be noted that our simulations follow the sub- 
mergence of the MF in the solid crust, as matter piles up 
on top of the crust, not the dynamical process of accretion 
outside. This would require MHD si mulations of the accr e- 
tion process in the exterior (see e.g. iBernal et all (2010)), 
which is out of the scope of this paper. Thus, we assume 
that any details about the dynamics of the accretion pro- 
cess is included in the phenomenological form of ra(0), that 
describes the rate at which matter reaches the top of the 
crust at each latitude. How the complex external dynam- 
ics leads to a particular form of m(0) should be the pur- 
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Figure 3. Magnetic field in the solid part of the crust after ac- 
creting 1O _4 M . Upper panel: model A with = 10 12 G and 
spherical accretion. Second, third, fourth panels: model A with 
B^ = 10 14 G, and (s), (p) and (e) accretion rate geometry, re- 
spectively. Fifth panel: model B with B® = 10 14 G and spherical 
accretion. 

pose of separate studies. We also note that at the relatively 
high temperatures (2 — 3 x 10 9 K) during the early accretion 
stage in our simulations, fast thermonuclear burning is likely 
to bring the composition close to the ground state, unlike 
what happens in accreting, old, cold NSs in binary systems. 
Therefore we assume that our crust is composed by matter 
in the ground state. 

We begin by discussing how our results depend on the 
value of M a for three representative cases: 
(s) spherically symmetric, ra(6) — k; 
(p) channeled onto a polar cap, rh(0) — ke~^°^° w ^ ; 

, 0-71-/2 ,2 

(e) channeled onto the equator, rh{0) — ke 1 °™ , 
where k is the normalization that fixes M a . Hereafter we 
discuss our results with 6 W = 0.4 rad. In Fig. [3] we compare 
the MF configurations in the crustal region immediately af- 
ter accretion stops, for models with M a — 10 -4 M©. We 
look first at the simple case of spherical accretion (upper 
two panels) for £° = 10 12 G and B® = 10 14 G for model 
A. During the accretion stage, field lines are advected to- 
ward the interior and screening currents are developed in 
the outer layers. The reduction in the surface field strength 
is compensated by the local amplification of Be at the sub- 
mergence depth, where it can reach values up to few 10 16 G. 
The geometry at the end of the submergence phase is not af- 
fected by the initial field strength (it simply scales with Bp), 
since the evolution is governed by the term in v a (Eq. Q), 
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Table 1. Properties for different accretion geometries and total 
accreted mass (M a ) for model A, with 5° = 10 14 G. The den- 
sity at which the maximum compression of lines is reached is pd- 
Bp and Bp stand for the values of B p at the end of the submer- 
gence stage and the maximum value reached during reemergence, 
respectively. 
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in Eq. (2]). However, anisotropic accretion flows cause the 
irregular submergence of the MF (middle and bottom pan- 
els). Strong screening currents also appear, but now they 
are localized in the regions where m{0) is larger. This has 
an important consequence: the external, large-scale, dipo- 
lar component is not reduced as much as in the spherically 
symmetric case. 

In Table [1] we summarize our results for model A, for 
different values of M a and varying the accretion geometry, as 
listed above. Here pd denotes the density at the submergence 
depth and Bp the final strength of the dipolar component 
of the MF (in all models B® = 10 14 G). We point out that 
the exact value of Bp is very sensitive to the boundary con- 
ditions, i.e. the magnetic flux exchange with the exterior. 
We match with a general vacuum solution, while a more 
consistent approach should consider again the MHD accre- 
tion process in the exterior. On the contrary, we found that 
our results for pd, Bp and reemergence timescales are not 
sensitive to the external boundary conditions. For spherical 
accretion (s), we found that if M a < 10 _5 M©, the submer- 
gence depth is shallow. We observe a sharp transition from 
M a < 10 _4 M©, with at most a factor 3 decrease in the sur- 
face strength, to M a > 1O _3 M0, characterized by a very 
deep submergence, well within the inner crust. 

Model B (fifth panel) shows a similar behaviour: the 
submergence depth is the same as model A. No differences 
in the submergence process arise from the geometry or in- 
tensity of the MF: the details of the M a — pd relation depend 
mainly on the particular structure and mass of the NS. As a 
matter of fact, in agreement with the induction and continu- 
ity equations, the enclosed mass between pd and the surface 
is M a , so that for M a > 1O _2 M0 the accreted mass has 
completely replaced the outer crust. 

The rotational behaviour of the star would also be rad- 
ically different, because the external dipole (third column) 
regulates the spindown of the star. For M a > 1O _3 M0, and 
spherical (or nearly isotropic) accumulation of matter on top 
of the crust, B p is largely reduced and therefore the spin- 
down rate of the star becomes unusually small. Although 
estimates of the MF based on timing properties suggest a 
very low field, the internal field could be orders of magni- 
tude larger, since it is simply screened by currents in the 
crust that will be dissipated on longer timescales. Note also 
that equatorial accretion (e) produces a similar reduction of 
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Figure 4. Reemergence of MF. Configurations at t = 5 kyr after 
the accretion episode, for the same cases shown in Fig. [3] 



the external dipolar field as for spherical accretion, but the 
submergence is rather anisotropic (deeper in the equator). 
On the contrary, accretion concentrated in the polar cap (p) 
barely modifies the MF strength or geometry. 

We have also found that adding an azimuthal compo- 
nent of velocity / in Eq. j4j) does not affect the sub- 
mergence process, driven solely by the radial inflow. Instead, 
its effect is the twisting of field lines and the creation of 
toroidal field, which in turn enhances the Hall activity on 
longer timescales, during the reemergence process. 

3.2 Reemergence of the magnetic field. 

After accretion stops, the reemergence process begins. It is 
mainly driven by Ohmic diffusion (screening currents are 
concentrated in the outer crust, where the resistivity is rela- 
tively high) , but also by the Hall term when the MF strength 
locally exceeds 2— 3x 10 14 G. In this latter case, there may be 
significant changes in the geometry of the MF, including the 
generation of higher order multipoles in both poloidal and 
toroidal compone nts, and the launching of whistler waves to - 
wards the poles (|Pons fe Geppertll2007l : IVigano et al]l2012h . 

In Fig. |4] we show the MF configuration at t = 5 kyr 
for the same five models as in Fig. [3] In the weak field case 
(top panel) the field simply reemerges by diffusion. It shows 
very little Hall activity, seen as generation of toroidal field 
(white/black spots) and deformation of poloidal field in a 
small equatorial region of the outer crust. In the other four 
cases (Bp — 10 14 G), the local generation of strong toroidal 
field and the creation of non-trivial structures are clearly 
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Figure 5. Radial profile at the equator of uos'Te for initial model 
A (left) and model B (right), with = 10 14 G, at three times: 
at the beginning (solid), after spherical accretion of 10 _4 Mq 
(dashed) and at t = 5 kyr, during reemergence (dot-dash). 



seen, but these small scale features are mainly localized in 
the inner crust, which makes difficult to predict possible 
observational consequences. At the same time, the global, 
large-scale field at the surface is being restored to a shape 
closer to the original geometry shown in Fig. [T] with addi- 
tional small scale multipoles near the surface. 

The importance of Hall activity is shown in Fig. [5] Here 
we plot the radial profile of the magnetization parameter at 
the equator, at different stages for model A (left panel) or 
model B (right panel), both with B° = 10 14 G. The solid 
line corresponds to the initial configuration, for which ujBT e 
is initially of order 1 for model A, while it is ten times larger 
for model B, due to the larger mean value of MF. After the 
submergence (dashed), in both models ujBT e reaches much 
larger values in the region where the field has been com- 
pressed. After 5 kyr, the diffusion of the field and the lower 
temperature make ujBT e > 1 in the whole crust. As ujBT e 
scales with P, the same models, but with Bp = 10 12 G, give 
a rescaling of ujBT e by two orders of magnitude, providing 
oJBTe < 1 always. Thus, for such initially low fields almost 
no Hall activity is expected. 

In the top panel of Fig. [6] we show the evolution of 
Bp as a function of time after a spherical fallback episode, 
with different M a , for model B (long dashes) or model A 
(other linestyles), together with the inferred values (or up- 
per limit) of observed CCOs. The reemergence timescale, on 
which the surface MF grows, depends basically on pd. For 
total accreted matter in the range of interest, there is an ini- 
tial delay of about 10 3 yr before appreciable reemergence is 
observed. The reemergence process takes between 10 3 and 
few 10 5 , and it is determined by the local conditions (in 
particular, the resistivity) where the screening currents are 
located. The MF will never reach the original strength, since 
some Joule dissipation is always expected. Generally speak- 
ing, Bp is restored to close to its initial value for the core- 
extended configurations (model A). Only for the extreme 
case of M a > 1O _2 M0, when the reemergence process lasts a 
long time, Bp is significantly reduced (see Table 1). In model 
B (long dashes), crustal currents are much stronger than in 
model A, so both Ohmic and Hall timescales are shorter. 
This implies that the reemerged field h as been more dissi- 
pated, in agreement with ID studies (|Geppert et al. 1999; 
1531201 ih . Furthermore, as the field is more compressed, the 
Hall activity is more intense (see dot-dash line in Fig. [SJ, 
and the transfer of magnetic energy from the dipolar com- 
ponent to higher multipoles (i.e., small scale structures) is 
faster. 
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Figure 6. Evolution of B p and P during the reemergence phase 
after spherical accretion, for Bp = 10 14 G, assuming Po = 0.1 s. 
Different cases are shown: model B, M a = 1O _4 M0 (long dashes); 
model A, M a = [1,5,10,20] x 10~ 4 M© (other lines from top to 
bottom, respectively). Observational data are shown with two 
romboids and an arrow (upper limit). Estimates and errors of the 
ages of the associated supernova remnants are shown 
and refs. within). 
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Figure 7. Comparison between surface temperature distributions 
of an anti-magnetar (solid), and two models of hidden magnetars 
(dashed and dot-dashed) after accreting M a = 10 _3 Mq. All three 
models are shown in an evolutionary phase when the external 
dipolar field is B p = 10 10 G. 

The reemergence stage has a strong imprint on timing 
properties. Assuming a small period at the end of the accre- 
tion stage, Po = 0.1 s, we follow the evolution of the timing 
properties (P and P) using the classical formula for an or- 
thogonal rotator: B p = (6.4 x 10 19 \/PPhh/Rlo) G, where 
Rio = R/10 km and I45 = //10 45 g cm 2 are given by the 
radius and moment of inertia of the star model. The bot- 
tom panel of Fig. [6] shows the evolution of P as a function 
of the real age for several values of M a in the sensitivity 
range [10 -4 — 1O _3 ]M0. Variations of a factor of a few in 
M a within this range result in an extremely low value of P 
during the first thousands of years of a NS life. In terms of 
the spin-down age, r c = P/2P, this means that it overes- 
timates the real age also by the same factor. Observational 
data reported for CCOs are consistent with M a > 1O _3 M . 
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A possible observable to distinguish between an "anti- 
magnet ax" and the hidden magnetic field scenario is the sur- 
face temperature distribution, T a (0), that determines the av- 
erage luminosity and light-curve in the X-ray band. Fig. 
shows T 8 (9) for two high field NSs during the reemergence 
phase after a 1O _3 M0 accretion episode and an antimag- 
netar. The dot-dashed line refers to model A (Fig. [T]) at 
t\ — 3 kyr, while the dashed line corresponds to a crust al- 
confined magnetic field model, with a larger toroidal com- 
ponent {B% = 10 13 G, B? = 10 15 G), at t 2 = 6 kyr. In both 
cases, the field has partially reemerged and has same value 
as for the antimagnetar model (solid lines), B p = 10 10 G. 
Model A has a slightly warmer surface than the antimagne- 
tar model, but both are almost isotropic. On the contrary, 
the crustal confined model shows a distinctive feature con- 
sisting in a lower average T s with hot polar caps. This is 
produced by the buried toroidal field that keeps the equa- 
torial surface region insulated from the warmer core. The 
anisotropy depends on the accreted mass, the natal mag- 
netic field and it varies with time depending on the location 
of currents during the reemergence process. 



4 FINAL REMARKS 

We have presented results from detailed simulations of one 
of the possible scenarios proposed to explain the low in- 
ferred external MFs of CCO. For the first time, a consis- 
tent, coupled magneto-thermal evolution in 2D and includ- 
ing the Hall term has been used to model the evolution of 
the magnetic field in both submergence and reemergence 
phase, under the effect of different accretion flow geometries. 
Our study c onfirms the qualitativ e results obtained by pre- 
vious works (|Geppert et al.|[l999h for ID simplified models, 
in what concerns the relation between submergence depth 
with total accreted mass. However, for large magnetic fields, 
the dynamics of the reemergence phase is more complex in 
our models, since we are not restricted to one particular 
component (dipole) and the interplay between toroidal and 
poloidal field is allowed. 

We have confirmed that the hidden magnetar model is a 
possible solution, but it requires that, immediately following 
the supernova explosion, 10 _3 M© are able to reach a large 
fraction of the the top of the crust. Under this circumstances, 
the accretion flow is able to push the field lines deep inside 
the crust. Other choices, with channeled accretion onto a 
polar cap or onto the equatorial region, have a more local 
effect that will barely modify the large-scale MF. If the NS is 
born as a typical pulsar (10 12 G), rather than a magnetar, a 
reduction of only 1-2 orders of magnitude of the ex ternal MF 
i s eas ier accomplished with less accreted matter. IChevalierl 
(|l989h suggested a value of M a ~ 0.1M© for SN 1987A, 
while M a < 1O _3 M is expected for a more typical type II 
SN. Interestingly, this estimate lies in our range of interest. 
Independent work modeling the MHD accretion process in 
the exterior of a NS surrounded by the dirty environment of 
the supernova explosion is needed to provide more realistic 
boundary conditions and to explain the required values of 
M a for a deep submergence of the field. 

We have also explored the subsequent evolution of the 
crustal MF during reemergence by Ohmic dissipation of 
screening currents. If the hidden MF is of the order of 10 14 G, 



non-linear dynamics driven by the Hall term are also ex- 
pected, with the generation of toroidal field and multipolar 
components. This activity is, however, localized in the in- 
ner crust, and it remains unclear if the surface field (and 
therefore timing properties) can reflect in some way this 
non-trivial structures. In any case, the reemergence process 
by diffusion has a clear imprint on the timing properties of 
the NS. The temporary low values of P may point to weak 
inferred MFs, but that is simply an artifact of the com- 
plex field geometry. The reemergence phase, for the explored 
range of M a , may last from 10 3 to 10 5 years, after which 
we expect that the field has almost completely reemerged, 
although it is partially dissipated. Another important ob- 
servational consequence of this scenario is the effect on the 
braking index, n. During reemergence, the field appears to 
be growing to an external observer, which means n < 3. 
Interestingly, this is in agreement with all reported values 
measured for young pulsars (< 10 4 yrs) ([Hobbs et al ] l2010l : 
lEspinoza et al.ll201lh , perhaps pointing to hypercritical ac- 
cretion in the aftermath of supernova explosions as a very 
common scenario. 

Last, the thermal evolution during the reemergence 
phase in the hidden magnetic field scenario can yield to 
totally different degrees of anisotropy in the surface tem- 
perature, depending on the initial MF and on the ac- 
creted mass. A hidden magnetar with a large, extended 
toroidal field buried in the crust could qualitatively ex- 
plain the surface temperatures anisotropics inferred in sev- 
eral CCOs: the large pulsed fraction observed i n Kes 79 
(jHalpern fe Gotthelj linol : IShabaltas feTal l2012h. the an- 
tipodal hot spots seen in Puppis A (|de Luca et al.|[201^ ). 
and the small emitting region of the b lackbody component s 
needed to fit the spectra of IE 1207 (|De Luca et al.ll2004l ). 
Instead, weaker (or less extended) hidden toroidal MF pro- 
vide low pulsed fractions. 

We conclude that the hidden magnetic field scenario 
for CCOs is not in contradiction with any observational 
data, and it would not be surprising that one of this ob- 
jects displays magnetar-like activity in the form of bursts or 
flares, as it has happened for the low-field magn etar SGR 
0418+5729 (|Rea et al.1 120I0I : iTurolla et afl l201lh . We also 
consider the possibility that not all the CCO candidates 
are of the same nature: large, magnetar-like hidden MFs 
are particulary suitable for those CCOs who are displaying 
temperature anisotropics. On the other hand, some of the 
standard, young rot at ion- powered pulsars (RPPs) could be 
in their final stage of reemergence, representing the connec- 
tion between CCOs and RPPs/magnetars population. 
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